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Abstract — This work presents joint iterative power allocation and in- 
terference suppression algoritlims for DS-CDMA networlis which employ 
multiple relays and the amplify and forward cooperation strategy. We 
propose a joint constrained optimization framework that considers the 
allocation of power levels across the relays subject to individual and global 
power constraints and the design of linear receivers for interference sup- 
pression. We derive constrained minimum mean-squared error (MMSE) 
expressions for the parameter vectors that determine the optimal power 
levels across the relays and the parameters of the linear receivers. In 
order to solve the proposed optimization problems efficiently, we develop 
recursive least squares (RLS) algorithms for adaptive joint iterative power 
allocation, and receiver and channel parameter estimation. Simulation 
results show that the proposed algorithms obtain significant gains in 
performance and capacity over existing schemes. 

I. Introduction 

Multiple collocated antennas enable the exploitation of the spatial 
diversity in wireless channels, mitigating the effects of fading and 
enhancing the perfoiTnance of wireless communications systems. Due 
to size and cost it is often impractical to equip mobile terminals with 
multiple antennas. However, spatial diversity gains can be obtained 
when terminals with single antennas establish a distributed antenna 
array through cooperation [r|-[3|. In a cooperative transmission 
system, terminals or users relay signals to each other in order to 
propagate redundant copies of the same signals to the destination 
user or terminal. To this end, the designer must employ a cooperation 
strategy such as amplify-and-forward (AF) |[3], decode-and-forward 
(DF) [3|, 141 and compress-and-forward (CF) 15]. 

Recent contributions in the area of cooperative and multihop 
communications have considered the problem of resource allocation 
0, |7|. However, there is very little work on resource allocation 
strategies in multiuser spread spectrum systems. In particular, prior 
work on cooperative multiuser spread spectrum DS-CDMA systems 
in interference channels has not received much attention and has 
focused on problems such as the impact of multiple access inter- 
ference (MAI) and intersymbol interference (ISI), the problem of 
partner selection |J4], JS] and the analysis of the bit error rate (BER) 
and outage performances J^. The problem of resource allocation 
and interference mitigation in cooperative multiuser spread spectrum 
systems has not been jointly dealt with so far. This problem is of 
paramount importance in cooperative wireless ad-hoc and sensor 
networks I.10J . 1.1 IJ that utilize spread spectrum systems. 

In this work, spread spectrum systems which employ multiple 
relays and the AF cooperation strategy are considered. The problem 
of joint resource allocation and interference mitigation in multiuser 
DS-CDMA with a general number of relays is addressed. In order 
to facilitate the receiver design, we adopt linear multiuser receivers 
114 1 which only require a training sequence and the timing. More 
sophisticated receiver techniques are also possible ( 14], 1 19 1. A joint 
constrained optimization framework that deals with the allocation of 
power levels among the relays subject to individual and global power 
constraints and the design of linear receivers is proposed. MMSE 
expressions that jointly determine the optimal power levels across 
the relays and the linear receivers are derived. Joint adaptive and 
iterative recursive least squares (RLS) algorithms are also developed 
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for solving the optimization problems, mitigating the effects of MAI 
and ISI, and allocating the power levels across the links. 

This paper is organized as follows. Section II describes a coop- 
erative DS-CDMA system model with multiple relays. Section III 
formulates the problem and the constrained linear MMSE design of 
the receiver and the power allocation. The proposed RLS algorithms 
for the estimation of the receive filter, the power allocation and the 
channels subject to a global and individual power constraints are 
developed in Sections IV and V, respectively. Section VI presents and 
discusses the simulations and Section VII draws the conclusions. 

II. Cooperative DS-CDMA System Model 

Consider the uplink of a synchronous DS-CDMA system commu- 
nicating over multipath channels with QPSK modulation, K users, 
A'^ chips per symbol and L as the maximum number of propagation 
paths for each link. The network is equipped with an AF protocol that 
allows communication in multiple hops using n,. relays in a repetitive 
fashion. We assume that the source node or terminal transmits data 
organized in packets with P symbols, there is enough training and 
control data to coordinate transmissions and cooperation, and the 
linear receivers at the relay and destination terminals are perfectly 
synchronized. The received signals are filtered by a matched filter, 
sampled at chip rate and organized into A4 x 1 vectors rsd[i\ and 
Vsri [i] which describe the signal received from the source to the 
destination and from the source to the relays, respectively, as follows 



rsd[n] = y^^a';a[n]Dkhsd,k[n]bk[i] 
fc=i 
+ 'nbd[n] +nsd[n], 

K 

-sTj [m] = 2J °-brj [m]Dkhsrj,k[m]bk[i] 
fc=i 

+ Vsr. ["l] + "sr, [r. 



(1) 



L2, 



imj, 
1 = 0,1, 



n = Hri + 1, m = n,.i + r j + 1 

where M — N + L — 1, risdli] and riarj [i] are zero mean complex 
Gaussian vectors with variance a^ generated at the receiver of the 
destination and the relays, and the vectors rj^^^ [i] and rj^^ . [i] represent 
the intersymbol interference (ISI). The M x L matrix Dk has the 
signature sequences of each user shifted down by one position at 
each column to form 



Dk 



dk{l) 

dk{N) 






dfc(l) 

dk{N) 



(2) 



where dk — [dfc(l), dfc(2), ..., dk{N)] stands for the signature 
sequence of user k, the L x 1 channel vectors from source to 
destination, source to relay, and relay to destination are hsd.k[n], 
hrd,k[n], hr-s,k[n], respectively. By collecting the data vectors in 



(TJ (including the links from relays to destination) into a (n,. + l)Mx 
1 received vector at the destination we get 



rsd[n\ 









ridr 



+ ri[i] + n[i] 
Rewriting the above signals in a compact form yields 

K 

r\i] = ^Cfc-Hfc[i]Bfc [i]aft[i] + 77(1] + n[{\, 



(3) 



where the {rir + 1) Af x (n,. + 1)L block diagonal matrix Ck contains 
shifted versions of Dk as shown by 



Ck — 



Dk 



D, 








Dk 



(4) 



The (rir + 1)L x {rir + 1) matrix Tikli] contains the channel 
gains of the links between the source and the destination, and 
the relays and the destination. The [ur + 1) x (n^ + 1) diag- 
onal matrix Bk[i] = diag(6fe[j] bj.^ [i] . . . &JI" [i]) contains the 
symbols transmitted from the source to the destination (bk[i]) and 
the Ur symbols transmitted from the relays to the destination 



(b'k 



I'l d r 



[i]) on the main diagonal, the (n,. + 1) x 1 vector 



ak[i\ = [cisdM i^ridi"^] • • • '^rn dM]'' of 'hc amplitudes of the 
links, the [rir + 1)M x 1 vector T7[i] with the ISI terms and 
{rir + l)Af X 1 vector n[i] with the noise. 

III. Problem Formulation and Proposed MMSE Design 

This section formulates the problem of MMSE joint power alloca- 
tion and interference suppression for a cooperative DS-CDMA net- 
work. Two constrained optimization problems are proposed in order 
to describe the joint power allocation and interference suppression 
problems subject to a global and individual power constraints. 

A. MMSE Design with Global Power Constraint 

The MMSE design of the power allocation of the links across the 
source, relay and destination terminals and interference suppression 
filters is presented using a global power constraint. Let us express 
the received vector in (|3) in a more convenient way for the proposed 
optimization. The {rir + l)Af x 1 received vector can be written as 



r[i] = CT'HT[i]BT[i]aT[i] + ri[i] + n[i]. 



(5) 



where the {rir + 1)M x K{nr + 1)L matrix Ct = [Ci C2 ... Ck] 
contains all the signatures, the K{nr + 1)L x K{nr + 1) block 
diagonal matrix "Wt [i] contains the channel gains of the all links, the 
K{nr + 1) X K{nr + 1) diagonal matrix Bk[i] = diag(fei[i] b'[^'''\i] 



.b\-%] 



bK[i\ fe^li] 



j;r„dr 



;]) contains the symbols trans- 



mitted from all the sources to the destination and from all the relays 
to the destination on the main diagonal, the K{nr + 1) x 1 vector 
otW = [aLW «ridM • • -"r^^dH, • • • , a^W a'^tdim] 
. . . a^ df'^]]"^ of the amplitudes of all the links. 

Consider an MMSE design of the receivers for the users rep- 
resented by a {rir + 1)M x K parameter matrix W\i\ — 
[wi[i], . . . , iVKli]] and for the computation of the K{nr + 1) x 1 
optimal power allocation vector aT,opt[i]- This problem is cast as 

[M^opt,aT,opt]=arg min E[\\b[i\ -W"[i\r\i]\f] 
subject to aT [i]aT[i] ~ Pt, 



where the K x 1 vector b[i] = [bi[i], 



, 6a'[J]]"^ represents the 



desired symbols. The MMSE expressions for the parameter matrix 
Wopt and vector ar.opt can be obtained by transforming the above 
constrained optimization problem into an unconstrained one with the 
method of Lagrange multipliers |20| which leads to 



£ = E[\\b[i\ - W" \i]{CTnT[i\BT[i\aT[i\ + r][i\ + n[i\)\ 
+ \{a?[i\aT[i]~PT), 



(V) 



Fixing arli], taking the gradient terms of the Lagrangian and 
equating them to zero yields 



Wopt = R ^Pc-H, 



(8) 



where the covariance matrix of the received vector is given by it = 

E[r[i\r"[i\] = CTnT[i\BT\i]aT[i\a^\i]Bf^[i\-H?[i\Cf^ + a^I 
and Pc^ = E[r\i]b"\i]] = E[CTnT\i]BT\i]aT[i\b''\i]] is the 
{ur + 1)M X K cross-correlation matrix. The matrices R and Pch 
depend on the power allocation vector arli]- The expression for 
arli] is obtained by fixing W[i], taking the gradient terms of the 
Lagrangian and equating them to zero which leads to 



flT.opt — {RaT + ^tI) Pa 



(9) 



where the K{nr + 1) x K{nr + 1) covariance matrix 
Rc^-r = E[Bi}[i]ni}[i]Ci}W[i]W"[i]CT'HT[i]BT[i]] and 
the vector p„^ = £'[Bf [i]-KTW"Cf W^bH] is a K{nr + l) x 1 
cross-correlation vector. The Lagrange multiplier At in the 
expression above plays the role of a regularization term and has 
to be determined numerically due to the difficulty of evaluating its 
expression. The expressions in (|8) and (|9) depend on each other 
and require the estimation of the channel matrix "HTli]- Thus, it 
is necessary to iterate (|8) and (|9) with initial values to obtain a 
solution and to estimate the channel. In addition, the network has to 
convey all the information necessary to compute the global power 
allocation including the filter Wopt. The expressions in ([8} and (|9) 
require matrix inversions with cubic complexity ( 0{{{nr + l)Af ) ) 
and 0{{K{nr + 1))^), should be iterated as they depend on each 
other and require channel estimation. 

B. MMSE Design with Individual Power Constraints 

Here, it is posed an optimization problem that considers the 
joint design of a linear receiver and the calculation of the optimal 
power levels across the relays subject to individual power constraints. 
Consider an MMSE approach for the design of the receive filter lUk [i] 
and vector ak[i\ for user k. This problem can be cast as 



[x(;fc,opt,afc,opt] = arg min E[\bk[i] 



-«;fWrH|^] 



subject to Ofe [i]afc[i] = P, 



(10) 



A,k, 



fc = l, 2, 



K. 



The expressions for the parameter vectors tyfc[i] and ak[i] can be 
obtained by transforming the above constrained optimization problem 
into an unconstrained one with the help of the method of Lagrange 
multipliers f20 1 which leads to 

K 

C = E[\bk[i]-w^[i]{Y,CiUi[i\Bi[i]ai[i\+r)[i\+n[i])\^] 
1=1 
+ \{a^[i\ak[i\-PA,k), k=l,2, ..., K. 

(11) 

Fixing Ofc [i], taking the gradient terms of the Lagrangian and equating 
them to zero yields 



Iffc.opt = R PcH^ 



1, 2, 



K, 



(12) 



where R = Ef=iCfc^feWBfcMafcHaf[J]Bf[J]?^f[i]Cf-fa2/ 
is the covariance matrix and p^^i = i5[&fe[i]r[i]] = Ck'H[i]ak[i] is 



the cross-correlation vector. The quantities R and Pc-h depend on 

ak[i]. By fixing ly^fj], the expression for ak[i] is given by 



Q-fcjopt 



= (i2„, +A/)-V , fc = l, 2, 



K, 



(13) 



where fl„, = J2LiBk[i]'Hkli\C"wk[i]w^[i\Ck'Hkli\Bk[i\ is 
the (rir + 1) X {rir + 1) covariance matrix and the (n^ + 1) x 1 cross- 
correlation vector is p„ = E[bk[i]B^[i]'Hk[i]"C"wk[i\]. The 
expressions in ill) and ( I13t depend on each other and should be iter- 
ated. They also require the estimation of the channel matrices Hk H] 
and need matrix inversions with cubic complexity ( 0{{{nr + l)M)' ) 
and 0{{nr + 1) ). In what follows, we will develop algorithms for 
computing a^^opt, 'Wk,opt and the channels 7ife[i] for k — 1, . . . ,K. 

IV. Proposed Joint Estimation Algorithms with a 
Global Power Constraint 

Here we present adaptive joint estimation algorithms to determine 
the parameters of the linear receiver, the power allocation and the 
chaimel with a global power constraint. 

A. Receiver and Power Allocation Parameter Estimation 

Let us now consider the following proposed least squares (LS) 
optimization problem 

i 

[Vt^W.aTM] =arg niin ^ a'-'||b[Z] - W^Hr[qf 

subject to a-j- [i]aT[i] = Pt, 

(14) 

where a is a forgetting factor. The goal is to develop a recursive 
cost-effective solution to J14t . To this end, we will resort to the 
theory of adaptive algorithms |20| and derive a constrained joint 
iterative recursive least squares (RLS) algorithm. This algorithm will 
compute W[i\ and aT[i] and will exchange information between 
the recursions. The part of the algorithm to compute 'W[i] uses 
*[i] = R[i] = J2Ui a'^'r[l]r"[l] and is given by 



k\i] = 



a ^^[i]r[i] 



(15) 



*[i] = a'^^[i - 1] - a''^k[i\r"[i\^[i - 1], (16) 

W\i]^W[i-l]+k[i]i"[i\, (17) 
where the a priori estimation error is given by 

^[i\ = b\i] - W" [i - l]r\i]. (18) 

The derivation for the recursion that estimates the power allocation 
ar [i] presents a difficulty related to the enforcement of the constraint 
and how to incorporate it into an efficient RLS algorithm. A modifi- 
cation is required in order to complete the derivation of the proposed 
RLS algorithm. This is because the problem in ( 1141 incorporates a 
Lagrange multiplier (A) to ensure the global power constraint, which 
is difficult to include in the matrix inversion lemma. 

Our approach is to obtain a recursive expression by relaxing 
the constraint and then ensure the constraint is incorporated via 
a subsequent normalization procedure. In order to develop the re- 
cursions for dT[i], we need to compute the inverse of Rar^li] = 
ELi-^tH'Ht WCTW[i]lV^[i]CT'HTWBTW- To this end, let us 
first define ^a^ = -Rqt [*] and employ the matrix inversion lemma 
1201 as follows: 

^ ^ a-^^c,T.[i\Bff[{\-HT[i\CffW[{\ 

1 + a-^w" [ilCT-HTlilBTlil^a^riilB^lilnTmcffWli] 

(19) 

^c,^[i\ = a-^^a.T.[i-l]-a-''ka.j,[i\w"[i]CTnT[ilBT[i\^c,j.li-'i-] 

(20) 



After some algebraic manipulations, we get 

arli] = &T[i] + fcarWCoT W (21) 

where the a priori estimation error for this recursion is 

Ca^M = bW - a?[i\B^[i]'HT\i]e?W[i\ (22) 

In order to ensure the individual power constraint afii[i]aT[i] = Pt, 
we apply the following rule 



otH ^ Pt aT[i][\/aT [i]aT[i 



(23) 



The algorithms for recursive computation of W [i] and ot [i] require 
estimates of the channel vector 'Ht [i] , which will also be developed 
in what follows. The complexity of the proposed algorithm is 

0{{{nr + l)Mf) for calculating W\i] and o{{K{nr + 1))^) for 
obtaining aT[i]- 

B. Channel Estimation with a Global Power Constraint 

We present a channel estimator that considers jointly all the K 
users and exploits the knowledge of the receive filter matrix V(^[i] 
and the power allocation vector aT[i\- Let us rewrite r[i] as 



r[i\ = CTBT[i\AT[i]hT\i] + r][i] + n[i]. 



(24) 



where Bt [i] is a Kijir + 1)1/ x K{nr + 1)L matrix with the symbols 
of all K users transmitted from the sources and the relays on the 
main diagonal, AT[i] is a K{nr + 1)L x K{nr + 1)L matrix with 
the amplitudes of all the links on the main diagonal and HtIi] is 
a Kijir + 1)1/ X 1 vector with the L channel gains from all K 
users, sources and relays. A channel estimator can be derived from 
the following optimization problem 



Ht \i] = are min >^ a 



CTBT[i]AT[i]hT[iW (25) 



The solution to the above optimization problem is given by 



where 



RhT.[i] = a Rh^[i-1]- 



hT[i\^RH.r\i\Ph^[^, (26) 



l + Q-l<[^]i^h^[^-l]^>T[^] 



Ph^ W = f^Ph^ [*-!]+ fT [i]r[i\ 



and 



(27) 
(28) 

(29) 



vrli] = CT-BT[i]AT[J] 
The channel matrix "HriJ] is determined with the recursion: 

UT[i\= Y. '^T\i]qT,j, (30) 

where qfj, ^ = [^_^^ 1 . . . ]. 

This algorithm jointly estimates the coefficients of the channels 
across all the links and for all users subject to a global power 
constraint. The complexity of the proposed RLS channel estimation 
algorithm is 0{(K{nr + 1)L)^)- 

V. Proposed Joint Estimation Algorithms with 
Individual Power Constraints 

In this section, we present adaptive algorithms to determine the 
parameters of the linear receiver, the power allocation and the 
channel subject to individual power constraints. Unlike the algorithms 
presented in the previous section, the algorithms proposed here are 
suitable for distributed resource allocation, detection and estimation. 



A. Receiver and Power Allocation Parameter Estimation 

We develop a recursive solution to the expressions in ( 1121 and 
U3b using time averages instead of the expected value. The proposed 
RLS algorithms will compute ty^fi] and ak[i] for each user k and 
will exchange information between the recursions. We fix dfc[i] and 
compute the inverse of R[i] using the matrix inversion lemma 11201 to 
obtain ii;fe[i]. Defining $[i] = R[i] then we can obtain the recursions 



The solution to the above optimization problem is 



m 



a ^^[i]r[i] 



1 + a^ir-'^[i]#[i]r[i] 



(31) 



*[i] = a"^*[i - 1] - a'^k[i]r"[i]^[i - 1] (32) 



Wk[i] =Wk[i-l\ + k[i\C[i] 
where the a priori estimation error is 

£,[i] = bk[i] -Wk[i- l]r[i]. 



(33) 



(34) 



The derivation for the recursion that estimates the power allocation 
follows a similar approach to the computation of liJfcfi]. In order to 
develop the recursions for ak[i]-, we need to compute the inverse of 
Rc.M = T:i=iB^[l\'H"[l]C^Wk[l]w^[l]Ck'hk[l]Bk[l\. To this 
end, let us first define ^a.,. = Ra^ [i] and proceed as follows: 



feo,. \i] 



Q-i*„Ji]Bf [j]?£,, [i]Cf Afc[J] 



1 + a-^w»[i\Ck'hk\i]Bk[i\^^^[i\B^\i\'k^[i\C^Wk\i] 

(35) 

(36) 

ak[i] = ak[i] + kc,,^ [i]^*^ [i] (37) 

where the a priori estimation error for this equation is 

Ca.M = bk[i\ - a^[{\B^[i\'H"[i\C^Wk[i\ (38) 

In order to ensure the individual power constraint a^[i]afe[i] — PA,k, 
we apply the rule 



fflfc I 



PA,k ak[i]{\l a" [i]ak[i] 



(39) 



The algorithms for recursive computation of ii!fc[i] and dfe[j] require 
estimates of the channel vector 'Kfc[i], which will also be devel- 
oped in what follows. The complexity of the proposed algorithm 

is 0(({nr + 1)A/)^) for calculating Wkii] and 0{{nr + 1)^) for 
obtaining ak.[i]. 

B. Channel Estimation with Individual Power Constraints 

We propose here an algorithm that estimates the channels for each 
user k across all links subject to individual power constraints and 
exploits the knowledge of the receiver filter Wk[i\ and the power 
allocation vector a,k[{\- Let us express the received vector in (3) as 



r[i\ =^CkBk[i\Ak[i]hk[i\ + r][i\ + n[i] 



(40) 



where Bk[{\ is a [rir + 1)1/ x (n,. + 1)1/ matrix with the symbols 
of user k transmitted from the sources and the relays on the main 
diagonal, A^fi] is a (n^ + 1)1/ x (n,. + 1)L matrix with the 
amplitudes of all the links on the main diagonal and hk[i] is a 
K{nr + 1)1/ X 1 vector with the L channel gains for user fc, and 
each link. A channel estimation algorithm can be developed solving 
the optimization problem 



/ife [j] = arg min y a ^\\r[l] - CkBk[i]Ak[i]hk[i]\\^ , 



h^[i] 






(41) 



hk W = Rh^ [i]p^,. W , for fc = 1, 2, ...,K (42) 



where 



Ry,^[i\=a-^Rf,l[i-l] 



^Rn^[i~l]vk[i\v'i;[mh^i-l] 
l + a-^v^[i]Rll[i-l]vk[i 



and 



Ph. W = aPfc. [*-!]+ Vk[i]r[i] 



Vh[i] = CkBh[i]Ak[i] 



The channel matrix "Hfefi] is obtained as follows: 

n,, + l 



(43) 
(44) 

(45) 



(46) 



where Qk.j = ^ 

This algorithm estimates the coefficients of the channels of each 
user k across all the links subject to individual power constraints. 
The complexity is 0{{{nr + 1)1/)^). 

VI. Simulations 

We evaluate the bit error ratio (BER) performance of the pro- 
posed joint power allocation and interference suppression (JPAIS) 
algorithms with global (GPC) and individual power (IPC) constraints 
and compare them with interference suppression schemes without 
cooperation (NCIS) and with cooperation (CIS) using an equal power 
allocation across the relays. We consider a stationary DS-CDMA 
network with randomly generated spreading codes with a processing 
gain A'^ = 16. The fading channels are generated considering 
a random power delay profile with gains taken from a complex 
Gaussian variable with unit variance and mean zero, i = 3 paths 
spaced by one chip, and are normalized for unit power The power 
constraint parameter PA,k is set for each user so that one can control 
the SNR (SNR = PA,k/o'^) and Pt = KPA,k, whereas it follows 
a log-normal distribution for with associated standard deviation of 3 
dB. We adopt the AF cooperative strategy with repetitions and all the 
relays and the destination terminal are equipped with linear MMSE 
receivers. We employ packets with 1500 QPSK symbols and average 
the curves over 1000 runs. The receivers have either full knowledge 
of the channel and the noise variance or are adaptive and estimate all 
the required coefficients and the channels using the proposed RLS 
algorithms with optimized parameters. For the adaptive receivers, we 
provide training sequences with 200 symbols placed at the preamble 
of the packets. After the training sequence, the adaptive receivers are 
switched to decision-directed mode. 

We first consider the proposed JPAIS method with the MMSE 
expressions of ((8) and (|9) using a global power constraint (JPAIS- 
GPC), and ( 112b and ( I13t with individual power constraints (JPAIS- 
IPC). We compare the proposed scheme with a non-cooperative 
approach (NCIS) and a cooperative scheme with equal power al- 
location (CIS) across the relays for n^ = 1,2 relays. The results 
shown in Fig. [T] illustrate the performance improvement achieved 
by the proposed JPAIS scheme and algorithms, which significantly 
outperform the CIS and the NCIS techniques. As the number of relays 
is increased so is the performance, reflecting the exploitation of the 
spatial diversity. In the scenario studied, the proposed JPAIS-IPC 
approach can accommodate up to 3 more users as compared to the 
CIS scheme and double the capacity as compared with the NCIS 
for the same performance. The proposed JPAIS-GPC is superior to 
the JPAIS-IPC and can accommodate up to 2 more users than the 
JPAI-GPS, while its complexity is higher. 



N-16, SNR - 12 dB, P=1500 




Number of users 



Fig. 1. BER versus SNR and K for the optimal linear MMSE detectors. 
Parameters: At = A = 0.025. 



N-1 6, SNR - 1 2 dB, K-6 users, P-1 500 




500 1000 

Number of received symbols 



Fig. 2. BER performance versus number of symbols. The curves for the 
adaptive schemes are in soUd lines, whereas those of the optimal MMSE 
schemes are in dotted lines. Parameters: Aj^ = A = 0.025 (for MMSE 
schemes), a = 0.998 (for adaptive schemes). 



The second experiment depicted in Fig. |2] shows the BER per- 
formance of the proposed adaptive algorithms (JPAIS) against the 
existing NCIS and CIS schemes with n,. = 2 relays. All techniques 
employ RLS algorithms for the estimation of the coefficients of the 
channel, the receiver filters and the power allocation for each user 
(JPAIS only). The complexity of the proposed algorithms is quadratic 
with the filter length of the receivers and the number of relays 
rir, whereas the optimal MMSE schemes require cubic complexity. 
From the results, we can verify that the proposed adaptive joint 
estimation algorithms converge to approximately the same level of 
the MMSE schemes, which have full channel and noise variance 
knowledge. Again, the proposed JPAIS-GPC is superior to the JPAIS- 
IPC but requires higher complexity and joint demodulation of signals, 
whereas the JPAIS-IPC renders itself to a distributed implementation. 

VII. Concluding remarks 

This work proposed joint iterative power allocation and inteifer- 
ence mitigation techniques for DS-CDMA networks with multiple re- 
lays and the AF cooperation strategy. A joint constrained optimization 
framework and algorithms that consider the allocation of power levels 
across the relays subject to global and individual power constraints 
and the design of linear receivers for interference suppression were 



proposed. The results of simulations showed that the proposed 
methods obtain significant gains in performance and capacity over 
existing non-cooperative systems and cooperative schemes. 
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